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We display some simple cosmological solutions of gravity theories with quadratic Ricci curvature 
terms added to the Einstein-Hilbert lagrangian which exhibit anisotropic inflation. The Hubble 
expansion rates are constant and unequal in three orthogonal directions. We describe the evolution of 
the simplest of these homogeneous and anisotropic cosmological models from its natural initial state 
and evaluate the deviations they will create from statistical isotropy in the fluctuations produced 
during a period of anisotropic inflation. The anisotropic inflation is not a late-time attractor in these 
models but the rate of approach to a final isotropic de Sitter state is slow and is conducive to the 
creation of observable anisotropic statistical effects in the microwave background. The statistical 
^\ • anisotropy would not be scale invariant and the level of statistical anisotropy will grow with scale. 



o 



> 

in 
o 



INTRODUCTION 



o 
cr, 

The observations of the cosmic microwave background (CMB) made by WMAP and ground-based detectors are in 
good general agreement with the expectations of a post-inflationary universe containing a residual quintessence field 
HQ. However, the data quality has led to a focus upon the presence of several unexpected features of the sky maps 
as well as a search for any evidence of non-gaussianity in their statistics In particular, there have been studies 
of the significance and possible explanations for a lower than expected power in the quadrupole signal Q , possible 
qq , alignments between low multipoles giving rise to a preferred direction, or 'axis', on the sky [5[ that can arise in some 
£f~) • anisotropic universes and an apparent asymmetry between the northern and southern ecliptic hemispheres Q. 
These deviations from statistical isotropy in the data are hard to assess definitively by means of a posteriori statistics 
and may be due to unnoticed biases in the foreground subtraction, but careful studies of this potential problem have 
yet to find evidence of such an effect at a level which can explain the observations Q. A detailed study of the 
evidence for homogeneous statistical anisotropy in the CMB has been conducted by Hanson and Lewis who use 
quadratic maximum-likelihood estimators to analyse Gaussian models with statistical anisotropy and realistic sources 
of instrumental noise. They find evidence for anisotropy in the power spectrum with a large angular dependence for 
the quadrupole, aligned close to the ecliptic plane. Since this is the plane in which the satellite moves there is a 
\ suspicion that this observed asymmetry may be associated with a systematic instrumental effect or a beam ellipticity 
■ - - i which is not corrected for in the CMB maps analysed. 

There have been several attempts to explore the possibility that small anisotropic features might be imprinted 
upon the primordial fluctuation spectrum by the process of inflation. In the standard general relativistic model 
of inflation driven by scalar fields which violate the strong-energy condition during a period of slow rolling, this 
would only be possible in the very earliest moments of inflation if the prior state was one of extreme expansion- and 
curvature-anisotropy because the dynamics rapidly approach those of the de Sitter metric in the presence of a positive 
effective cosmological constant. This situation is studied in ref [Toj . but only for the simplest form of expansion 
anisotropy ignoring the effects of collisionless particles 11 1 and curvature anisotropies [13, Il3j|. which both make the 



residual effects larger. It is also possible to induce local expansion anisotropies by the presence of super- horizon scale 
inhomogeneities, and this scenario is explored by Gao [l4| ■ If we change the underlying inflationary model by adding 
a vector field (l5l-[l7j. a Chern-Simons or Kalb-Ramond field [18], or some quadratic curvature corrections to the 
gravitational lagrangian then the situation can change. Effective stresses are created which (unlike the scalar 
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inflaton field in general relativity) can violate the dominant or the weak-energy conditions as well as the strong- 
energy condition. Under these circumstances the de Sitter metric may no longer be an attractor during accelerated 
expansion and anisotropic inflationary behaviour might even be possible. Anisotropic inflation is only possible in 
general relativity with a positive cosmological constant and matter that obeys the strong energy condition if the 
3-curvature is positive. There is a known exact S 1 x S 2 Kantowski-Sachs universe of this type [20] which expands 
forever, with exponential expansion of its scale factors in two directions and a constant scale factor along the third. 
This model has been used in a study of inflation [2lj . however, its anisotropic behaviour is unstable within this class 
of Kantowski-Sachs models and to the addition of matter fields 22]. 

In this paper we are interested in the situation where a quadratic Ricci term, R I1V R} IV , is added to the Einstein- 
Hilbert term. We have already shown [3] that in the resulting theory with positive cosmological constant it is possible 
to find simple exact spatially homogeneous cosmological models of Bianchi type I which inflate anisotropically. Here 
we explore these solutions in more detail and show that the same cosmological models also possess an exact de 
Sitter attractor solution. However, the evolution towards the de Sitter attractor from an earlier state of anisotropic 
inflation is unusually slow and there is ample opportunity for anisotropic statistical effects to be imprinted upon the 
fluctuation spectrum created by the accelerated expansion. These anisotropics will be larger at earlier times in the 
inflationary phase and therefore will imprint greater anisotropic effects on larger scales than on smaller ones. This 
model is mathematically simple, with Euclidean space sections, and provides a useful testing ground for computing 
more detailed effects of anisotropic inflation on the scalar and tensor irregularity spectra. It allows us to determine 
the sense and relative magnitudes of explicit expansion anisotropics during a period of volume inflation in which the 
orthogonal scale factors expand by different exponential time factors. 

Usually inflation is considered to be driven by an isotropic scalar field; here we will consider a simple cosmologial 
constant as an example of such an effect. The models can easily be extended to include isotropic scalar fields of this 
kind, however, this may require additional fine-tuning for the scalar field to come to dominate the expansion dynamics 
at the right time. This is also the case for f(R) effects through the standard conformal transformation [25|. The 
effects from the quadratic Ricci term, on the other hand, are different as they allow for anisotropic degrees of freedom 
- indeed, a quadratic Ricci term is the main source of the anisotropic inflation considered in this paper. 



II. A BIANCHI TYPE I MODEL 



Our staring point is the equations of motion using a dynamical systems approach. For the Bianchi type I (and II) 
models, these are given in [241 ] . We will consider the quadratic theory where the Einstein-Hilbert action is modified: 



S G = J- / d 4 x^/\g~\ (R + aR 2 + PR^R^ - 2A) 
2^ Jm 



= 2aR [R^ - -R 9fiV ) + (2a + 0) (g^D - V M V„) R 

+/3D ( Rfj, v - \Rg^ ) + 2/3 ( R^ vp - - A g^R ap ) R ap , (2) 



i M 

Upon variation, the equations of motion are the modified Einstein equations: 

G>„ + + Ag^ = 0, (I) 
where = R pv — {\/2)Rg^ u is the regular Einstein tensor, and 

I 

4 

^Rg^ + 2/3 (r^ vp - J; 

and □ = V M V M . We note that the GR-limit can be obtained by letting (a,/3) — > (0,0). Now consider the spatially 
homogeneous Bianchi type I metrics. We can always write their metric line-elements as 

ds 2 = -At 2 + 5 ab uj a uj b , 

where uj a is a triad of one- forms which, for the Bianchi type I model, can be written as: u) a = e a i (t)dx l . 

Defining the time-like hypersurface-orthogonal vector u = d/dt, we can define the Hubble scalar, H, and the shear, 
(Tab, as follows: 

H = °ab — U( a - b ) - HSab- 

We will also restrict attention to cosmological models where the shear is diagonal, so we can write 

a a b = diag(-2cr + , a + + \/3a-, a + - v3<7_). 
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We define the dimensionless expansion-normalised variables by scaling out appropriate powers of H 

R 1 & 

B = 7Z 1 mrr9 > X 



{'3a + (3)H 2 ' 3a + /3 7 

H H A 



ff 2 ' ^ if 3 ' 3# 2 ' 
E± = ^> E±i = §h S ±2 = gj.. (3) 

Note the presence of time derivatives of the variables Q 2 and £±2; this reflects the 4* ft -order time derivatives in the 
field equations of the quadratic theory [13]; X is a constant. We also introduce the dynamical time variable r by 

At 

and note that since H = a/ a, the dynamical time is related to the length scale as a = ioe T , where i o > is a constant. 
We also assume that the cosmological constant is positive: f2 A > 0. 

The full set of dynamical equations for the Bianchi type I model (and type II model) are given in [24j . In [24| it was 
shown that there exists a peculiar set of solutions for this theory, namely a set of anisotropically inflating solutions. A 
stability analysis of these solutions was also performed. Here, we intend to study these solutions further; in particular, 
we will investigate their possible influence on any statistical relic anisotropy at the end of inflation. 

For later reference, let us also give the equations of motion. We will, for simplicity, consider the invariant subspace 
when (£+, £_) oc E_i) oc (£+2, £-2)- This enables us to write: 

(E+, £_) = £(cos<^, sin0), (S+i, £-1) = £1 (cos 0, sin0), 

(E +2 ,E_ 2 ) = E 2 (cos0,sin0), cj>' = 0. (4) 
The equations of motion are now: 

B' = -2QB, (5) 

n' A = -2Qn A , (6) 

Q' = -2Q 2 + Q 2 , (7) 
Q' 2 = -3(Q + 2)Q 2 -^(Q + 2)Q-^b(i + Z 2 -n A + ^Q 

-^(l + 2 X )S 4 -i(8 + x)S 1 -(4-x)SS 1 

-i(4- X )(3£ 2 + 2£S 2 + 2Q£ 2 ), (8) 

£' = -QE + Ei, (9) 
E; = -2QEX + E2, (10) 

E' 2 = -3(Q + 2)E 2 + ^i [B-(llx-8) + 4Q(l-x)+4E 2 (l + 2 X )] 
X 

+ - [W + (4 - x )(6 + O2 + 7Q) + 4(1 + 2 X )(3E 2 + 2EEi)] (11) 

These equations are subject to the ('Friedmann-like') constraint: 

= B(l - n A - E 2 ) + 12Q - 2Q 2 + 4Q 2 - (4 - x)(3 + 2Q)E 2 

-6(1 + 2 X )E 4 - X (E 2 - 2EE 2 ) + 4(2 + X )EEx. (12) 

The parameter Q is related to the usual deceleration parameter q via 

q = -(l + Q), 

and the variable B measures how greatly the quadratic part of the lagrangian dominates over the general-relativistic 
Einstcin-Hilbcrt term R — 2A. In particular, the larger the value of B, the "closer" we are to GR. The B — case 
corresponds to a purely quadratic lagrangian theory whose equations of motion reduce to $ M „ = 0. 
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A. The Inflating solutions 

We will now focus on two (sets of) solutions, namely, the de Sitter solution, and the anisotropically inflating type 
I solutions. 

a. The de Sitter solution, dS: The de Sitter solution is characterised by the critical points where 



E±i - e ±2 = N = o, n A = i, B^O. 



Its stability is assured if [24 



B >0^(3« + /3 )>0, fl + 2(4- X ) <0=> l + 2A(4a ±R <0 . 

X P 

b. Anisotropically-inflating type I universes, A(I): For certain values of x an d B, there are also exact solutions 
that describe anisotropic inflationary solutions of Bianchi type I [24[ : 

(E+,E_) = E(cos^,sin0), E 2 = E 2 = 2(4-x)+5 



4(2* +1) 



E±i = E ±2 = N = 0. 



There are two classes of such solutions, depending on the values of B and f2\- Here we will concentrate on the 
following: 

B = constant, il\ = 



8(2* +1) 



So long as x an d B take values for which Eq > 0, these solutions exist. Moreover, the solution is unstable, but for 
B > contains only one unstable mode. 

The metrics corresponding to the case where B ^ can be written 



d.s^ = -dr + e 



2 , 26* 



- 4CT +*d:r 2 + e 2(- + + V3--)' d J, 2 + e 2(^+-V3a_)t d ^2 



(13) 



2 1 + 8A( Q + /?) ,2^2^ 1 + 2A(4o + 0) 



We note that for these solutions to exist we need both these squares to be positive. Furthermore, we require er 2 + 
a\ + b 2 /2 = A/3 > 0. 

III. ANISOTROPIC INFLATION 

Based on the solutions above and their stability, there are values of the parameters where A{I) exists, and the de 
Sitter solution is stable. In fact, the connection is deeper than this. Assuming B > 0, and A > 0, then if % > 4, 
A(I) is connected to dS via a bifurcation at B — 2(x — 4). For B < 2(x — 4) the de Sitter solution is unstable (1 
unstable mode). As the value of B increases, dS bifurcates at B — 2(x — 4) creating the equilibrium points A(I). 
For B > 2(x — 4), A(I) has acquired the unstable mode and made dS stable. Hence, the only unstable mode of A{I) 
is actually connected to the de Sitter solution and these inflationary solutions should therefore be connected via a 
heteroclinic orbit (at least close to the bifurcation point). It is this fact we shall exploit here because this opens up 
the possibility for that the universe approaches the point A(I) and so starts to inflate anisotropically. This state is 
almost stable but it does have an unstable mode. This unstable mode will therefore eventually drive the evolution 
away from anisotropic inflation, towards isotropic inflation (represented by the stable point dS). 

In Figure [T] the evolution of universes with three different initial values are plotted. The universes start to inflate 
around r = 0, for which Q w 0, however, as can be seen, the shear is non-zero; hence, the universe starts inflating 
anisotropically. During a transient period the universe can stay arbitrary close to the anisotropically inflating point 
A(I); however, as we can see, eventually the universe will move towards isotropic inflation. For the universes displayed, 
we see that during the time from r = to r = 60 (which corresponds to 60 e-folds), we have Q ~ and therefore they 
are inflating during the entire period, even if the universes are undergoing a transition between anisotropic inflation 
and isotropic inflation. 

Interestingly, the smaller the shear is at r = 0, the longer it takes to enter the state of isotropic inflation. The 
reason for this can be seen from the unstable mode of A(I). From the eigenvalues derived in (24|, we see that the 
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FIG. 1: Inflation: The evolution of Q and E with time, r, for universes with three different initial values. Here, we set x — 6. 
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unstable mode goes as oc exp(Air), where Ai = (3/2) (-^/l + 8E 2 , ~ !)• Consequently, the smaller the shear, the more 
time it would take for the universe to enter a state of isotropic inflation. Furthermore, since the universe inflates 
during the transition, Q and the value of B ss constant. Therefore, close to the de Sitter isotropic point the shear 
will decay as £ oc exp(A27~) where A2 can be estimated to be 

For small E i this reduces to A2 ~ — 4(2x + 1)Sq/(3%). This implies that the shear during the anisotropic inflation, 
E , is imprinted in the decay rate of the shear on the approach to isotropic de Sitter state and the smaller the shear 
So, the slower the approach towards the isotropic de Sitter inflation. 

A consequence of this is that at the end of 60 e- folds (r = 60), the universe may still have considerable amount of 
shear in spite of the fact that the universe has inflated. However, as we can see from Figure [TJ this shear can also be 
arbitrary small. 

In Figure [2] the numerics have been started away from the anisotropic inflationary solutions. With a bit of fine 
tuning we can see that the solutions experience a transient period during which the universe inflates anisotropically. 
This shows that there is a set of non-zero measure of initial values that experience anisotropic inflation. 

A. Consequences for the CMB 

In ref. [l6l |. the imprints of a preferred direction on the CMB was studied for an anisotropic model in which 
rotational invariance is broken by the presence of a vector picking out a preferred direction with unit vector n. If 
parity is preserved (k — > —k) then the leading order of the anisotropic power spectrum has a quadrupole form, with 

P(k) = P(k)[l + g(k)(n ■ k) 2 + higher order], 

where P{k) is the isotropic part of the power spectrum and g(k) measures the power of the statistical anisotropy. If 
the anisotropic contributions are scale independent then g(k) will be a constant. However, this need not always be 
the case, as in the study of anisotropies induced by super-horizon inhomogeneities . 

The model considered in ref (l6j was an axisymmetric anisotropically inflationary model. Our model allows for 
axisymmetry in the special case £_ = (the general case has no rotational symmetry whatsoever). Their model 
actually corresponds to the exact equilibrium point, but in our model this is unstable and evolves toward the isotropic 
de Sitter state. In their eq.(39) they define the anisotropy parameter en , which in our notation is just the dimensionlcss 
shear distortion (in the special case £_ = 0): 

g(k) oc£ ff E Hf± = 2S, 

(recall that the dynamical time r is related to the scale factor a through a — aoe r ). In our model g(k) will evolve 
and so cause the perturbation spectrum to depend on the time when each perturbation scale left the horizon. Hence, 
the statistical anisotropy will not be scale invariant. In our models the anisotropy is larger at the beginning of the 
inflation which means that long wavelengths will be more anisotropic than the short wavelengths. 

B. Other Bianchi models 

The model considered here is a very special Bianchi model. However, the existence of anisotropic inflation seems 
to be a widespread feature of quadratic theories of gravity. In [23| we showed the existence of anisotropic inflation 
for Bianchi type II models. Among the most general Bianchi models are the Bianchi type VI ^ and VII/j models and 
anisotropic inflation is also present there. The following metrics can be shown to be solutions to $ M „ = in cq. (fTJ): 

ds 2 = -dt 2 + dx 2 + e 2r{t+x) [e 2a (Ady + Bdzf + e- 2a (Cdy + Adz) 2 ] . 
Here, r and a are constants. For the various Bianchi types, the functions A, B and C are as follows: 



1. Type Vll h : A = cos[u(x + 1)], B = -C = sm{uj(x + t)]. 
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FIG. 2: The evolution of the dimensionless expansion-normalised shear, E, with time, r, for three universes starting away from 
the inflating solutions. They experience a transient period during which the universes inflate anisotropically. The evolution 
displayed requires some fine-tuning but the plot illustrates that there is a set of non-zero measure that gives the desired 
behaviour. 



2. Type Vl h : A = cosh[w(a; + 1)], B = C = smh[uj{x + t)]. 

3. Type IV: A = 1, B = lo(x + t),C = 0. 

There are no analogues of these solutions in general relativity because they would violate the weak energy condition. 
Interestingly, these are also so-called plane-wave spacetimes. We note that the 3-space volume expands exponentially, 
V cx exp(2ri) in terms of the comoving proper time, t\ hence, these are indeed inflationary solutions. However, as we 
see that one orthogonal direction is actually fixed and does not expand at all. 

The existence of such solutions for these Bianchi types indicates that anisotropic inflation is a more general feature 
of quadratic theories than previous thought and is present in even the most general class of anisotropic universes, 
which includes types Vlh and Vllh- 



IV. CONCLUSIONS 



Simple Bianchi type I universes in gravity theories in which the Emstein-Hilbert lagrangian is augmented by 
the addition of terms quadratic in the scalar curvature (R 2 ) and Ricci invariant {R liV Rf J,v ) display evolution that 
commences from a near-isotropic singularity. It mimics the behaviour of a radiation-dominated Friedmann universe 
and asymptotes towards a de Sitter late-time attractor. However, the evolution spends a long time evolving 
slowly through a period of anisotropic inflation during which the 3-volume expands exponentially in comoving proper 
time and the three directional scale factors increase at different rates. We display exact solutions which display this 
transitional anisotropic inflationary behaviour and show that it does not arise in the general relativity limit when the 
higher-order curvature terms vanish from the lagrangian. With a small amount of fine tuning of the pre-inflationary 
evolution this behaviour can leave a distinctive imprint on the spectrum of inhomogeneities created by any period of 
inflation defined by the exponential increase in the 3-volume. Using the characterisation of the leading statistically 
anisotropic contribution to the power spectrum of inhomogeneities that was introduced by Ackerman, Carroll and 
Wise [16( | , we determine the amplitude of the statistical anisotropy in the spectrum and show that it will not be scale 
invariant because the anisotropic effects increase with scale. 

There have been a number of studies of the possible sources of statistical anisotropy in the power spectrum of the 
microwave background. The most detailed study @ finds significant evidence for anisotropy but no proof that it is 
primordial in origin. Our analysis identifies a broad class of higher-order gravity theories in which there are solutions 
that differ significantly from those in general relativity. In particular, ever-expanding vacuum solutions with a positive 
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cosmological constant do not approach the de Sitter solution. The higher-order curvature terms contribute effective 
stress terms which violate the energy conditions that are needed for the cosmic no hair theorems of general relativistic 
cosmology to hold. Hence, they allow new types of exact solution to exist in which different directions accelerate 
at different rates. We have shown that solutions of this anisotropic inflationary type also exist in some of the most 
general Bianchi type universe and are not confined to the type I case we have used for simplicity. In the general 
Bianchi I universe we have studied, there is attraction to an asymptotic de Sitter solution, as in general relativity, 
but the evolution spends a large number of e-folds of expansion in the neighbourhood of an anisotropic inflationary 
solution. The existence of such behaviour near the Planck scale and the persistence of anomalies in the sky maps 
derived so far from observations of the microwave sky suggest that there is a possibility the two are related. 
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